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We investigate the relationship between the extensor electromyogram (EMG) and tremor time series 
in physiological hand tremor by cross-spectral analysis. Special attention is directed to the phase 
spectrum and the effects of observational noise. We calculate the theoretical phase spectrum for a 
second order linear stochastic process and compare the results to measured tremor data recorded from 
subjects who did not show a synchronized EMG activity in the corresponding extensor muscle. The 
results show that physiological tremor is well described by the proposed model and that the measured 
EMG represents a Newtonian force by which the muscle acts on the hand. 



1 Introduction 

Time series of hand tremor and the related muscle ac- 
tivities of the flexor and extensor muscles are obtained 
by measuring the acceleration of the hand (denoted here 
by ACC) and the surface electromyogram (denoted here 
by EMG). The ACC data of physiological tremor have 
been described as a linear stochastic process driven by 
uncorrelated firing motoneurons (Stiles and Randall 1967; 
Randall 1973; Rietz and Stiles 1974; Elble and Koller 
1990; Gantert et al. 1992; Timmer et al. 1993). The 
description of physiological tremor by a linear model is 
reasonable because linear approximations hold due to its 
small amplitude. These linear stochastic processes and 
their spectral and cross-spectral properties were studied 
exhaustively (Bloomfield 1976; Brockwell and Davis 1987; 
Priestley 1989). Usually, they are denoted by autoregres- 
sive processes, since actual values are given by a linear 
combination of past values plus a driving noise. In terms 
of physics, these processes are linear damped oscillators 
driven by noise. 

In the context of linear stochastic processes the rela- 
tion between two processes can be analyzed by investigat- 
ing phase and modulus, i.e. coherency, of the normalized 
cross-spectrum. Applications of cross-spectral analysis to 
EMG and ACC data of physiological tremor are reported 
in (Fox and Randall 1970; Pashda and Stein 1973; Elble 
and Randall 1976; Stiles 1983; Iaizzo and Pozos 1992). 
Up to now, the coherency and the phase spectrum were 
investigated only at a single frequency. In particular, the 
phase was always interpreted as a time delay between the 
two processes. 

The interpretation of the phase spectrum as a whole is 
difficult. For example, as will be shown below, the phase 
spectrum between EMG and ACC time series depends 



only on the mechanical properties of the hand and does 
not allow to draw conclusions on the dynamics of the driv- 
ing force, i.e. the EMG. In general, the phase spectrum 
can only be interpreted under quite strong assumptions 
about the interrelation of the processes. We discuss those 
cases relevant for the EMG - ACC relationship. Finally, 
we compare the spectra predicted from the model with 
those estimated from measured data. 

The paper is organized as follows: In the next section 
we briefly describe the data. In Section || we introduce 
the mathematical background for this and a companion 
paper (Timmer et al. 1998). Section [| gives theoretical 
and empirical results for physiological tremors showing a 
flat EMG power spectrum, resulting from unsynchronized 
muscle activity. EMG power spectra exhibiting a synchro- 
nization and the possible role of reflexes are discussed in 
a companion paper (Timmer et al. 1998). 

2 The data 

The data were recorded from normal subjects. The 
recording technique is described in detail elsewhere 
(Deuschl et al. 1991). Briefly, the time series of the hand 
tremor (ACC) were measured by a light-weight piezo- 
resistive accelerometer. The sampling rate is 300 Hz. 
The outstretched hand is supported at the wrist. We 
recorded three data sets for each subject, the first with 
the hand unloaded, the second with a 500gr load and the 
third with a lOOOgr load. The weights were fixed on the 
belly of the outstretched hand. External elements as the 
amplifiers and the piezoresistive sensors produce additive 
white observational noise in each recorded time series, un- 
correlated to the measured dynamical process itself. The 
variance of the observational noise can be estimated from 
the high frequency part of the power spectrum where the 
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contribution of the tremor oscillation can be neglected. 
This noise contributes up to 10% of the variance of the 
recorded data and has a significant effect on the estimated 
coherency and phase spectra as will be shown in Section 

I 

The EMG time series (EMG) were measured by sur- 
face electrodes fixed over the belly of the extensor carpi 
ulnaris muscle and the flexor carpi ulnaris muscle. These 
data represent broad band noise. The information about 
a possible synchronization of the muscle activity is en- 
coded in a modulation of this noise. The data were high 
pass filtered (cut-off frequency 80 Hz) in order to remove 
movement artifacts, rectified in order to obtain time series 
reflecting the muscle activity (Journee 1983) and then low 
pass filtered (cut-off frequency 150 Hz) to avoid aliasing. 
Finally, the signals were digitized and fed into a computer 
for off-line analysis. 

Like the tremor time series the EMG time series are 
contaminated with additive observational noise. Its vari- 
ance cannot be estimated analogously to that of the ACC 
data, since uncorrelated EMG activity also shows a flat 
power spectrum at higher frequencies indistinguishable 
from that of the observational noise. 

Time series from 58 subjects who showed no signifi- 
cant EMG synchronization, i.e. a flat spectrum, were ex- 
amined. The statistical decision of consistency with a 
flat spectrum was performed by means of a Kolmogorov- 
Smirnov test at the level of confidence p=0.05 (Timmer 
et al. 1996). A representative example of such a time se- 
ries is shown in Fig. [l]. Time series from 19 subjects with 
enhanced physiological tremor who showed a significant 
EMG synchronization are analyzed in the companion pa- 
per. In each time series the mean was subtracted and all 
series were scaled to variance one. 
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Figure 1: Acceleration of the hand (a) and rectified EMG (b) 
of physiological tremor. 



3 Mathematical methods 

In this section we introduce the mathematical methods 
that will be used in following sections and the companion 
paper to analyze the simulated and the measured data. 
Firstly, we briefly summarize the time and frequency do- 
main properties of linear stochastic processes before dis- 
cussing the cross-spectral estimation and interpretation. 
Special attention will be payed to the effects of observa- 
tional noise which is always present in the data of physio- 
logical tremor and its EMG and renders the interpretation 
of the mathematical results more difficult. 

3.1 Linear stochastic processes 

An example of a linear stochastic process is the autore- 
gressive (AR) process of order p : 

v 

x(t) =5^o,-aj(*-t) + e(f) , (1) 
i=i 

where e(f) denotes uncorrelated Gaussian noise with vari- 
ance a 2 . For ease of notation we set the sampling interval 
to unity for the theoretical discussions. Such a process 
can be interpreted in terms of physics as a combination of 
relaxators and damped oscillators (Honerkamp 1993). For 
example, an AR process of order 2 with appropriate pa- 
rameters a\ and ai describes from a physical standpoint 
a linear, damped oscillator with characteristic period T 
or frequency lj — 1/T, and relaxation time r. T and r 
are related to the parameters a\ and a,2 by: 

01 = 2cob0^ exp(-l/r) (2) 
a 2 = -exp(— 2/r) . (3) 

AR processes can be generalized to the autoregressive 
moving average (ARMA) processes by including past driv- 
ing noise terms in the dynamics. A more substantial gen- 
eralization is the linear state space model (Honerkamp 
1993). It allows one to model explicitly the observational 
noise 77(f) that covers the dynamical variable x(t) which 
is mapped to the observation by C and contributes to the 
measured z(t): 

x(t) = Ax(t-l) + e(t) (4) 
z(t) = Cx(t)+ri(t) . (5) 

This model has been applied successfully to physiological 
tremor time series (Gantert et al. 1992). If the observa- 
tional noise is not modeled explicitly, e.g. by applying an 
ARMA model, the characteristic times will be underesti- 
mated and statistical tests to decide on the model order 
will fail to detect the correct order (Konig and Timmer 
1997). This might explain the high model order reported 
earlier (Randall 1973; Miao and Sakamoto 1995). 
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3.2 Spectral properties of linear stochas- 
tic processes 

The power spectrum S(w) of a mean zero and unit vari- 
ance process x(t) is defined as the Fourier transform of 
the autocorrelation function ACF(r) : 

ACF(r) = < x(t) x(t -t)> (6) 

SH = ^^ACF(r) e — ,L,e[-n,7r] (7) 

Z7T * — ' 

T 

with "< >" denoting expectations. The estimation of the 
power spectrum is usually based on the Fourier transform 
F(uj) and the periodogram Per(w) of the data : 
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X(w) = -^=^2x(t) exp(-iwt) 
Per(w) = 
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and is evaluated at the frequencies: 
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The expectation of the periodogram is the power spec- 
trum but the periodogram is not a consistent estimator 
for the power spectrum since the standard deviation of 
this \\ distributed random variable is equal to its mean 
and does not decrease with increasing number of data 
(Brockwell and Davis 1987; Priestley 1989): 



PerH ~ 1 -${lu) X 1 



(11) 



In order to estimate the power spectrum, the periodogram 
has to be convolved by a window function W(J) of width 
2/i+l : 
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It is also possible to estimate the power spectrum by aver- 
aging the periodograms of segments of the data or by fit- 
ting an AR process to the data and calculate the spectrum 
of the fitted process. General aspects of spectral estima- 
tion as well as confidence intervals are given in (Brockwell 
and Davis 1987; Priestley 1989). Special aspects concern- 
ing spectral estimation for tremor time series are discussed 
in (Timmer et al. 1996). 

For linear processes the power spectrum can be calcu- 
lated analytically. In the case of an AR process of order 
2 it is given by: 



S(u) = 
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2ir |1 — ai exp(— iuj) — ai exp(— 2iw) 



(13) 



Expressed in terms of T and r, the power spectrum shows 
for | cos(27r/T)| cosh(l/r) < 1 a peak at the frequency: 



Wpeafe = arccos (cos(27r/T) cosh(l/r)) 



(14) 



Therefore, for small r, the peak of the power spectrum is 
not located at the frequency 2n/T. The width of the peak 
is proportional to 1/r. If the driving force is characterized 
by some nontrivial power spectrum Sd r ive(^>) instead of 
the constant spectrum of uncorrelated white noise ( filf ) 
changes to: 
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|1 — ai exp(— ibS) — 02 exp(— 2iuj)\^ 



(15) 



3.3 Cross-spectral analysis 



Analogously to the univariate quantities introduced in the 
previous section the cross-spectrum CS(w) of two zero 
mean and unit variance time series x(t) and y(t) is de- 
fined as the Fourier transform of cross-correlation func- 
tion CCF(t) : 

CCF(r) = <x{t)y{t-r)> (16) 

CS(w) = — VCCF(r)expHa;r) (17) 

Z7T L — ' 

r 

= <X{uj)Y*{uj)> . 

Here * denotes complex conjugation. The coherency spec- 
trum Coh(cj) is defined as the modulus of the normalized 
cross-spectrum CS(cj) : 
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(18) 
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and the phase spectrum $(w) by the representation: 

GS(w) = |CS(w)| eaq>(i#(«)) . (19) 

It can be shown that Coh(u;) equals one whenever y(t) 
is a linear function of x(t). This holds especially for the 
coherency between an AR process and its driving noise 
e(t) (Brockwell and Davis 1987; Priestley 1989). Thus, 
the coherency can be interpreted as a measure of linear 
predictability. The interpretation of the phase spectrum is 
more difficult. For the following cases the phase spectrum 
can be calculated analytically : 

- If the process y(t) is a time delayed version of process 
x(t), i.e. y(t) — x(t—At), the phase spectrum is given 
by a straight line with its slope determined by At : 



(20) 



If y(t) is the derivative of x(t), i.e. y(t) = x(t) a 
constant phase spectrum of — ^ results. 



= -tt/2 



(21) 



In the case of an AR process of order 2 (AR[2]), the 
phase spectrum between the driving noise e(t) and 
the resulting process is given by : 



$(w) = arctan 



a\ sin lu + a-2, sin 2uj 
1 — ai cos to — a 2 cos 2lo t 



(22) 
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It is important to note that the phase spectrum does 
not change if the driving noise is not a Gaussian white 
noise process. Because of the linearity of the sys- 
tem, eq. ( p2| ) holds whenever the driving process 
shows a broad band power spectrum. It might even 
be chaotic. If the relaxation time t is not smaller 
than the period T, the phase relation $di scr (cj) be- 
tween an (timc-discrctc) AR[2] process and its driv- 
ing noise ( p2[ ) is in good approximation related to 
the well known phase relation $ co „ t (a;) for a (time- 
continuous) differential equation of a linear, driven 
damped oscillator by: 



(23) 



Eq. (|23|) shows that there is a substantial difference 
between a discrete- and continuous-time treatment 
of the data, since modeling continuous-time data by 
discrete time models yields to an spurious time delay 
of one unit of the sampling period. Although, in the 
case of tremor data, the natural approach would be 
the continuous- time version, we choose the discrete- 
time version because the simulation studies become 
much easier and the mentioned effect is easily cor- 
rected for. 




Figure 2: Power spectrum of an AR[2] process (a). Phase 
spectrum between the AR process and its driving noise (b). 
The solid line gives the result according to the dashed 

line the result for the continuous time case taking ( fjq j into 
account. 

Fig. ||a displays the power spectrum of an AR[2] pro- 
cess and Fig. ^|b the phase spectrum between the AR 
process and its driving noise. The solid line gives the 
result according to (p2|), the dashed line shows the re- 



sult for the continuous time case taking (|23| ) into account 



Fig. g demonstrates that an interpretation of the phase 
spectrum $(w) for a single frequency is only possible un- 
der strong assumptions about the relation between the 
processes. In particular, the interpretation of the phase 
spectrum at a single frequency ujq, e.g. the frequency of 
maximum coherency, as a time delay by At = $>(u>q)/wo 
may be erroneous. This situation is similar to the in- 
terpretation of power spectra where in general a certain 
amount of power at a certain frequency may not be inter- 
preted as an oscillator of this variance. The whole phase 
spectrum, on the other hand, can provide substantial in- 
formations on the relation between the processes if the 
empirical phase spectrum fits to one of the theoretical 
phase spectra ( p0| , ^l| , p2] ) or combinations of them. 

The cross-spectrum CS(cj) is estimated analogously to 
(|l2|). The critical value s for the null hypothesis of zero 
coherency for a significance level a is given by: 



1 - a"- 2 



(24) 



where v is determined from the window function W(j) 
by: 

v = —, . (25) 

Confidence intervals for the coherency are given in Bloom- 
field (1976). Besides the simple case where y(t) and x(t) 
are indeed uncorrelated, at least the following reasons can 
result in a coherency unequal one: 

- A nonlinear relationship between x{t) and y(t) 

- Additional influences on y{t) apart from x(t) 

- Estimation bias due to misalignment (Hannan and 
Thomson 1971) 

- Observational noise 

If y(t) is a linear function of x(t) but the measurement 
of y(t) is covered by white observational noise of variance 
tj 2 oh the coherency is given by (Brockwell and Davis 1987) 
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Thus, the coherency is a function of the frequency depen- 
dent signal to noise ratio. For the general case of obser- 
vational noise on both processes the coherency is given 
by: 



Coh(w) = a 1 



{S x +al){S v + aD 



(27) 



where the argument ui was suppressed on the right hand 
for ease of notation and er^ and a 2 denote the constant 
power spectra of the observational noise. Fig.|| illustrates 
( p7f ) for different signal to noise ratios of both processes. 
This might partially explain the findings of Stiles (1980) 
and Lenz et al. (1988), who report a correlation between 
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peak power and coherency at the peak frequency as an 
effect of observational noise. If we assume a constant 
amount of observational noise the peak power is correlated 
with the signal to noise ratio. 
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Figure 3: Coherency between linear related processes in the 
case of observational noise on both processes for different sig- 
nal to noise ratios (SNR). The abscissa displays the SNR for 
one of the signals. The different curves parameterize the SNR 
for the other process. 



Eq. (|27]) is of particular interest since the variance of 
the estimator for the phase spectrum <£>(w) depends on 
the coherency (Priestley 1989): 
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where v is given by (|25|). Eq. 28 holds if the coherency 
is significantly larger than zero. For a coherency towards 
zero, the distribution of the estimated phase approaches 
the uniform distribution in [— tt,tt]. Therefore, the phase 
spectrum cannot be estimated reliably in the case of small 
coherency. 

Using (28), theoretical phase spectra like ( ptj| ) or (f22| ) 
can be fitted to estimated phase spectra by a maximum 
likelihood procedure. We used the Levenberg-Marquardt 
algorithm (Press et al. 1992). This algorithm provides 
confidence limits for the estimated parameters, that are 
asymptotically valid. The asymptotic results hold in the 
finite case when the parameter estimates are Gaussian. 
In order to test whether this condition applies in our case 
we performed a Monte Carlo simulation for an AR[2] pro- 
cess under conditions analogous to those observed in the 
empirical data. The variance of the driving noise, i.e. the 
unsynchronized EMG, was set to unity, the frequency of 
the AR process was 10 Hz, the relaxation time 0.1 s, 
i.e. one period. For a sampling frequency of 300 Hz, ac- 
cording to (@,[|) the parameters of the AR[2] process are 
ai = 1.8922 and a-i = —0.9355. Gaussian observational 
noise was added to both processes to obtain a signal to 
noise ratio of 10:1. Fig. |] shows scatter plots the esti- 
mated frequency, relaxation time r and delay Ai for 500 



realizations of the process. The Gaussianity of the esti- 
mates is clearly visible. Furthermore the estimates are 
uncorrelated. This is plausible since the period T deter- 
mines the frequency at which the phase spectrum is vastly 
varying, whereas the relaxation time t is related to the 
steepness of the phase spectrum at the that frequency. 
The estimated delay time corresponds to the sampling 
period of 0.0033 s due to (p3|). The variances and the co- 
variance of zero are consistent with the results from the 
Levenberg-Marquardt algorithm. Note that the relative 
error in r is much larger than that in T . The goodness- 
of-fit is judged by the \ 2 statistic (Press et al. 1992). 
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Figure 4-' Scatter plots of estimated frequency and relaxation 
time t (a) and delay At and relaxation time r (b) from the 
Monte Carlo study. 

The peak frequencies of the power spectra were esti- 
mated by the frequency of maximum power. A bootstrap 
method to obtain confidence regions for the true peak fre- 
quency will be described elsewhere (Timmer et al. 1997). 
Briefly, many periodograms are simulated from the esti- 
mated power spectra by the relation (11) and the power 
spectrum is re-estimated. The quantiles of distribution of 
the peak frequencies from the re-estimated power spectra 
yield a confidence region for the estimated peak frequency. 

If two linear processes with autocorrelation functions 
ACFi(i) and ACF 2 (t) are independent, the estimated 
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cross-correlation function CCF(t) is Gaussian distributed 
as : 

(5CF(r)~A/ , ^0,JV- 1 ^AC3Fi(i)ACF2(i)J • (29) 

Therefore, independence is difficult to infer from the cross- 
correlation function since the confidence interval depends 
on the autocorrelation functions of both processes that 
are in general not known. Only if one of the processes is 
white noise, the 95% confidence interval for a zero cross- 
correlation function is given by ±1.96JV — . Furthermore, 
the estimated cross-correlation function is not uncorre- 
cted for different lags. The covariance is given by : 

(cCF(ti)CCF(t 2 )) = l^ACF 1 (i)ACF 2 (t + r 2 -r 1 ). 

(30) 

Again, the autocorrelation functions of both processes en- 
ter the equation (Brockwell and Davis 1987). If, for ex- 
ample, one process is white noise and the other is oscillat- 
ing the cross-correlation function will show an oscillating 
behavior, even if the processes are independent. If the 
processes are not independent, ( |29| , |30| ) contain also the 
true cross-correlation function (Bartlett 1978), rendering 
the interpretation even more difficult. 

In the following section we apply the methods intro- 
duced above to measured data of physiological tremor 
without synchronization in the EMG. Whether and how 
cross-spectral analysis can contribute to decide about the 
origin of synchronized EMG in the case of (enhanced) 
physiological tremor is discussed in a companion paper 
(Timmer et al. 1998). In both cases the flexor EMG 
appeared to have a negligible contribution to the ACC 
data, i.e. the coherency spectrum is most often consistent 
with the hypothesis that the processes are uncorrelated. 
Whenever there was a significant coherency it was invoked 
by cross talk from the extensor EMG. The amount of cross 
talk can be estimated from the discontinuity at lag zero 
of the cross-correlation function because of its instanta- 
neous effect. The dominant contribution of the extensor 
is plausible since it is the anti gravity muscle. Thus, only 
the extensor EMG is considered in the analysis. 

4 Results 

It was frequently observed that a tremor appears even 
without synchronization in the EMG. This was inter- 
preted as a resonance phenomenon and described by an 
AR process (Stiles and Randall 1967, Randall 1973, Gan- 
tert et al. 1992). Fig. [| shows the results of the spectral 
and cross-spectral analysis for the data displayed in Fig. |l|. 
Fig. ||a shows the corresponding spectra, Fig. ||b the co- 
herency spectrum, Fig. ||c the phase spectrum and Fig. ||d 
the cross-correlation function estimated as described in 
Section 3. The straight line in Fig. ||b represents the 5% 
significance level for the hypothesis of zero coherency. The 



dashed line in Fig. [|d gives the 5% significance level for 
the hypothesis of zero cross-correlation assuming that at 
least one of the processes is white noise according to (pS|). 
The phase spectrum is shown 2tt periodically for a range 
of ±37r. The confidence regions of 2tr are only given for 
the central curve. 

The coherency spectrum seems to exhibits two peaks at 
approximately 8 and 12 Hz. Taking the confidence regions 
for the true coherency into account which are not shown 
for sake of clarity reveals that these peaks are not signif- 
icant, but represent a single peak in the region 7 to 14 
Hz. The fact that the coherency spectrum shows its max- 
imum values in the region of the peak of the ACC power 
spectrum, can be explained by (|2^) since EMG and ACC 
data are contaminated with noise. Compared to Fig. ^ 
the phase spectrum of the data is shifted by ir. Since 
we measured the acceleration instead of the position this 
results from applying (|2l|) twice to the ACC data. For fre- 
quencies below 3 Hz the small coherency and, therefore, 
the large errors of the estimated phase disable its inter- 
pretation. Note that the cross-correlation function shows 
a periodic structure also for negative time lags. Although 
they are statistically not significant, one could speculate 
whether they give evidence for some kind of reflex feed- 
back. 

In the frame of linear stochastic processes (with- 
out reflex feedback) the different spectra and the cross- 
correlation function should be determined by the follow- 
ing six parameters. 

• The characteristic times T and r, determining the 
parameters a\ and a 2 of the AR process modeling the 
mechanic properties of the musculosceletal system. 

• A possible time delay At between the EMG an ACC 
data. 

• The variance varEMG of the white noise e(t) model- 
ing the asynchronous EMG activity. 

• The variances var i, s . acc and var b s , emg of the ob- 
servational noises rait). 

Denoting the EMG by y(t), the movement of the hand 
by x(t) and the measured values by the subscript m the 
model reads : 

y(t) = e(t), e(t) ~ N{Q,var E M G ) (31) 
x(t) = at x(t - 1) + a 2 x(t - 2) + y(t - At) (32) 
y m (t) = y(t)+ m (t) (33) 

rji (t) ~ Af(0, var obs . emg) 
x m (t) = x(t)+r] 2 (t) (34) 
r) 2 (t) ~ Af(0, var obs . AC c) 

By (glg^^gg^) we fitted the parameter to the 
data. First, we fitted the phase spectrum without taking 
a possible time delay into account. This resulted in an 
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Figure 5: Results for a physiological tremor without EMG 
synchronization, (a): power spectra (EMG: dashed line, ACC: 
solid line), (b): coherency spectrum, the straight line repre- 
sents the 5% significance level for the hypothesis of zero co- 
herency, (c): phase spectrum with 95% confidence intervals, 
(d): cross-correlation function, the dashed lines display the 5 
% significance level for zero cross-correlation assuming that at 
least one of the processes in consistent with white noise. Con- 
fidence intervals for the power spectra and the coherency are 
not displayed for reasons of clarity. 



inappropriate fit. Only the inclusion of a time delay ac- 
cording to (|23| ) into the model gave a fit consistent with 
the data. Note, that this time delay of one sample unit 
does not reflect a time delay between the processes under 
investigation. It results from using a time-discrete model 
to describe an originally time-continuous process as dis- 
cussed in Section 3.3. A realization of the fitted model 
and the estimated spectra are displayed in Fig. [| Taking 
into account the errors of all estimated quantities, it shows 
good quantitative agreement with the empirical results of 
Fig. ||. The decreasing coherency for the high frequencies 
due to the frequency dependent signal to noise ratio and 
the resulting errors in the phase spectrum are well repro- 
duced. For the low frequencies the coherency of the model 
seems be larger than that of the data. This phenomenon 
was confirmed in many data sets. The discrepancy be- 
tween the data simulated by the model and the measured 
data results from the contribution of the heart beat to 
physiological tremor (Elble and Koller 1977). This addi- 
tional influence on the ACC apart from the EMG is not 
captured by the model. Thus, the coherency of the mea- 
sured data is reduced more than expected from the model 
that only considers the effect of observational noise. 

From a comparison of Fig. ||d and ^d it can be con- 
cluded that the small oscillations of the cross-correlation 
for negative lags give no evidence for a reflex feedback. 
These oscillations appear because the cross-correlation es- 
timates are not uncorrelated as discussed in Section 3.3. 

Assuming the validity of the AR[2] model to describe 
the physiological tremor, one can compare the peak fre- 
quency estimated from the power spectrum (14) with that 
estimated from the phase spectrum according to (p|,p|]22|) . 
Taking the errors of the estimates into account both val- 
ues are consistent. 

We received similar results in 70% of the investigated 
series. In the other 30% of the cases an interpretation of 
the phase spectrum was not possible because of the poor 
coherency and, therefore, large errors in the phase spec- 
trum. As discussed in Section 3 this might be simply the 
result of a smaller signal-to-noise ratio in the EMG and/or 
ACC in these cases due to a very low tremor amplitude. 
An estimation of the power in the ACC and EMG spectra 
supported this hypothesis. 



5 Summary 

We investigated the relation between muscle activity and 
physiological tremor in the case of unsynchronized EMG 
activity by cross-spectral analysis with special respect to 
the phase spectrum and the effects of observational noise. 
Such an analysis is not a straightforward task since one 
can not a priori expect a one to one relationship between 
one muscle and a specific mechanical measure like force, 
movement or acceleration due to the redundancy of the 
muscle system. Furthermore, the analysis is handicapped 
by the small tremor amplitude. 
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We found that this type of physiological tremor can 
be regarded as an AR process driven by the uncorrelated 
EMG activity without involving any reflex mechanisms. 
We showed that the phase spectrum between EMG and 
ACC can not be interpreted at a single frequency in terms 
of a delay. The phase spectrum depends on the mechani- 
cal properties of the hand, i.e. the driven part of the sys- 
tem, but not on the characteristics of the driving force. 
The behavior of the coherency spectrum can be explained 
as an effect of the frequency dependent signal-to-noise ra- 
tio. In addition, for low frequencies, the effect of the heart 
beat on the tremor further reduces the coherency between 
the EMG and the ACC. 

Autoregressive processes of order 2 are derived from 
stochastic differential equations where the noise repre- 
sents a force in a Newtonian sense, i.e. causing an acceler- 
ation. The conformity of the theoretical phase spectrum 
assuming such a process with the empirical data shows 
that in the case of this small amplitude hand tremor the 
measured EMG represents a Newtonian force by which 
the muscle acts on the hand. 
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